NBI-HE-93-06 
DFTT 1/93 
April 1993 
hep-lat / yymmnn 



COLOURED POLYMERS 



Igor PesandoQ 

The Niels Bohr Institute 
Blegdamsvej 17, DK-2100 Copenhagen 
Denmark Q 



Abstract 

We show that non-oriented coloured polymers (self-avoiding walks 
with different types of links) are in the same universality class of the 
ordinary self-avoiding walks, while the oriented coloured are not. 
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In this letter we would consider the critical behaviour of coloured poly- 
mers, both oriented and unoriented. A coloured polymer is a polymer with 
different kinds of monomers that have different binding energies among them; 
these coloured polymers will be described, following De Gennes's approach 
(Ijl]]), by a self-avoiding walk with different colours on links. 

We are interested in considering such models because: 

• 1) the unoriented coloured polymers model is not obtainable with a 
generalization of De Gennes model but only generalizing the Hil- 
horst one And before the analytic continuation down to n = 0, 
the Hilhorst model and the generalized one belong to different univer- 
sality classes: the model, we present, has nC components (hereafter C 
indicate the number of colours) with global symmetry Snc® ^2 while 
Hilhorst model has n components with Sn ® TL^ symmetry. Hence it 
is not, a priori, completely clear whether the analytic continuation to 
n = will bring the two model in the same universality class as one 
could expect naively. 

• 2) as we will explain below, the model we consider gives meaning to 
negative statistical weights. 

• 3) they seem to be more realistic that De Gennes' model since 
in nature there are polymers that are not composed by one kind of 
monomer only. 

• 4) the coloured oriented polymers model can be in a different class of 
universality with respect to the usual polymers. 

The model, we consider for unoriented polymers, is defined by the follow- 
ing partition function: 

c ^ 

Zn = Trexp{K ^ Sai ■ AijSbj Wab) (1) 

i,jS:G a, 6=1 

where K is the temperature of the model, Aij is the adjacency matrix of the 
lattice (or more generally the graph ) Q Sai = (S^i) is a n-components vec- 

-•^At least these two generalized models are not equivalent off criticality as we will see. 
^ Sn is the symmetric group of n elements. 

■^We use i,j,k... to indicate a site of the lattice f/, a, /3, 7 ... to denote the components 
of the vector S 
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tor defined on the site i with colour a, {w)ah = Wab is the matrix of statistical 
weights that determines the probability of the transition from colour a to 6; 
we point out that Wab = Wba- The trace is defined as follows 

Tr{*)= I n \{dSa.^Y.^{S:,"-^C) n 5(5fj . (2) 

■' ieg a=l "'-^ a,a l3^a,a 

This expression means that only one component of the vector {Su . . . Sci) = 
(BaZi Sai is different from zero ant its value can be ±a/ nC; hence the average 
on a site i is done on 2nC configurations characterized geometrically by the 
fact that they are the centers of the faces of a hypercube of side in 
nC dimensions. 

Following Hilhorst we get immediately that in the limit n ^ 

<S:X> = ^'^^^ab 

<S:y...SZ> = Vn>2 (3) 
It is also easy to find that 

limZ„ = l (4) 

n— >0 

and that the two points function yields at n = 

<S:,S^^>=6^^Y.Na^■,bAL) (5) 

L 

where Nai-bj{L) is the number of self-avoiding walks that have ends in i and 
j where there are colours a and h respectively; explicitly we have 

Na^,bAL) = Ni,AL) {w'^)ab (6) 

Would we have used the De Gennes approach, i.e. if we would have used 
as density p{Si) oc Wa^iS^i ~ ^) the definition of the trace, we still had 
found but we had not found because, differently from the one colour 
De Gennes model, here not only the self-avoiding walks would contribute 
to @ but also C tolerant walks (i.e. walks that can self intersect at most 
C times): this happens because we have different colours, and, for instance, 
the trace Tri^S^g^S^^Si^^S^fj) oc (1 — Sab) is different from zero using De Gennes 
approach and it allows for configurations like: 
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nevertheless the critical behaviour of C tolerant coloured walks is the same 
of the coloured self avoiding walks exactly as it happens to usual SAW and 
k-tolerant walks(0). 

In order to prove this affirmation, we notice that any symmetric real 
N X N matrix can be diagonalized by an SO{N) rotation, in particular we 
can choose a i? G SO{C) such that 

w = R^ A R (7) 

where A = diag{Xi . . . Ac) is the matrix of eigenvalues of w, that can be 
either positive or null or negative: with this meaning this model gives an 
interpretation to negative weights also ( see also eq. (Qj ). Redefining 

S'ia = R^ia (8) 

we can immediately rewrite as 

Z„ = Tt exp(if E >^aS',^ ■ A A) (9) 

i,j&g a=l 

Before we can perform a gaussian transform on this expression we have to 
take care of the antiferromagnetic couplings using the standard way: firstly 
we divide the lattice Q (supposed not frustrated) in two sublattices Q'^^'^ and 
Q^'^ in such a way each site of the first (second) sublattice can be reached 
by an even number of steps from the sublattice itself and by an odd number 
from the second (first) sublattice, then we fiip the sign of all the spins in the 
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first sublattice getting in this way a ferromagnetic coupling^. If we perform 
the gaussian transform we get 

Zn = cst / n n d$ai e"^ ^^^ee $a.-A^-$a,+i:. v;,,(v^0O (10) 
where in the Hilhorst case the effective potential is 



1 nC 

i + ^ElAai(C)^ + ^EA^(C)^ + ... (11) 



a, a 



while in the De Gennes case it turns out to be 



DG 



n i+iEiA.i(c)^+^As(B«y^) + 



(12) 

At the critical dimension we can consider only the relevant operators and the 
effective actions become^ 



S, 



eff H 



ba ■ (-A + ml)$a + Ua{<flf + n VaJ2 ' 



)4 

aa 



SeffDG = f d''xJ2[$a-{-A + ml)$a + ga{&] (13) 

that in the limit n — ^ yield the same result because the cubic anisotropy 
coefficient is proportional to n. Moreover at the critical point only the fields 
with zero mass propagate: these have the maximum absolute value of the 
eigenvalues Xm = max{| Aa |} since oc — || |. We will denote by 
J-" = {a\l < a < C s.t. I Aa I = Xm} the set of the indices of the critical fields 
and with JF^ the subsets J-'^ = {a\l < a < C s.t. Aa = ±Aa./}. 

^ The difference between AF and F coupiing can be found in the fact that for AF 
coupling the magnetic field is not uniform but has opposite values on the two sublattices 

g^+^ and g^-'> 

^Notice the coupling constants related to the null eigenvalues of A are obviously zero. 
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Let us now reconsider the asymptotic behaviour of the 2-point functions 
that can be now rewritten as 

< ^ai^bj > = X! ^paRpb < ' S^j ' >~ ^ RpaRpb < <Ppi ' (ppj ' > 

^pgjF+ RpaRpb + ( — J2pe:F- RpaRpb 



6-^ (E, 



(14) 



where we used the well known fact (0) that both F and AF share the same 
critical exponent r]. From what follows we immediately deduce that the 
coloured polymers have the same critical exponents of usual SAW. There is 
another way of getting the same result and it consists on considering the 
generating function of all the coloured polymers having a fixed end: 

Xab = H < ^ai^bj >= RpaRpb < Spj ' >= ^ RpaRpb ^"'^ ^(^) ^p)^ 

j P j PL 

(15) 

the first singularity happens when KXm = where z is defined by the 
asymptotic behaviour N{L) ~l>>i z^L^^'^, and hence we get 

x^f = E < s:,s^j >~ '^"''E ( E Rpc^Rpb + (-)^ E Rp^Rp] {^kxm)'^l''-' 

j L \peJ^+ p(^T~ j 

(16) 

We now show that if the matrix w is symmetric and not block diagonal 
then all the fields become critical to the same "temperature" , i.e. Va, h I]pgjF RpaRpb 7^ 
0. This is easily demonstrated if we note that a symmetric matrix w that 
is not block diagonal is irreducible^, the vectors = [Rpa) at fixed -p are 
eigenvectors of w and we apply the Perron-Frobenius theoremQfor the non- 
negative matrices (0]): in particular it follows that the sets are made 
just of one element. 

^ In the subsequent text a vector with positive (non-negative) coordinates will be called 
a positive (non-negative) vector, and similarly for a matrix. 

A matrix A is called reducible if there is a permutation matrix P such that P^^AP 
X Y 
Z 

irreducible. 

^ The part of the Perron-Frobenius theorem, we need, states that: An irreducible 
non-negative matrix A always has a positive eigenvalue r that is a simple root of the 
characteristic polynomial. 

The modulus of any other eigenvalue does not exceed r. 



fx Y 

is of the form _ 'z 1 1 where X and Z are square matrices. Otherwise, A is called 
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One could still wonder whether a model for oriented self-avoiding walks, 
that could admit an asymmetric weight matrix w, could exhibit also different 
critical behaviour: the answer is yes. An asymmetric weight matrix can 
be generated if we allow for different binding energies for each of the two 
(different) ends of a monomer: for instance, the weight matrix for a one- 
colour oriented model can be written as 

^=f^++ ^+-') (17) 

where ± distinguish between the two possible orientations of the monomer 

within the oriented polymer, and we used the fact that w = as it is 

shown by the following picture where ei , 62 are the two different endpoints 
of the monomer 

ei 62 Ci 62 62 61 62 61 



10++ — w 

The model for oriented coloured SAW is defined by the partition function 
(a generalization of (P) obtained through the approach of (0)) 

/ C7 \ 

Zoriented Coloured SAW = Tr exp I K ^ '^{daiAijPaj + IJaiA-ijlJaj) j (18) 

where aai, Pai are commuting fields and rjai, f]ai are grassman fields. The non 
vanishing traces are, by definition: 

Trl = 1 

TrUaPb = Trrjafib = Wab (19) 
Using these definitions it is easy to show that 

Zoriented Coloured SAW = 1 
< 0,aiPbj >= ^^Nai-bj[oriented){L) (20) 



To the "maximal" eigenvalue r there corresponds a positive eigenvector. 
Moreover, if A has h eigenvalues of modulus r, then these numbers are all distinct and 
are roots of the equation \^ — = 0. 
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where Nai-bj (oriented) (L) is analogous to the formula (||). If we perform a 
gaussian transformation on (18), we find 



Oriented Coloured SAW 



. C 

/ n n ^^ai ^'^'^i ~ iKWab{uJai(pbi + i^aiAi] 



a=l ieG 



(21) 



where ujaj and (f)aj are bosonic fields defined in j e ^ with "oriented" colour 
a and similarly for the grassmann fields ipaj and i/jaj- It can be easily shown 
that near the criticality 



Wa 



< (f)ciUJdj > WdbOcYl N, 



ai\bj (oriented) 



(L) 



(22) 



In the model defined by (^) we can use both asymmetric and symmetric w, 
but differently from the symmetric case where by means of a field redefinition 
we can always put the action in a diagonal form, in the general asymmetric 
case we cannot diagonalize the quadratic part of the action: this happens 
because an asymmetric real matrix can at most put in the Jordan canonical 
Ji 

R^^ where the generic Jordan block is of 



form, i.e. w = R 



the form X. 




. If we perform the naive continuum 



limit of (|2l|) and we throw away the most massive fields, we get 



Seff = J d^xiY^ (uJa{-A5ab + Mab)(j)b + tfja^-^^ab + Mab)i^b)-{ Gabii^a4>b + Ipai'i 
a,b \ a,b 

(23) 

where M and G are two upper triangulated invertible matrices with all the 
diagonal element equal R. 



^Notice that this action is stable under the symmetries 

1 
P 



uj puj , 
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Since a positive matrix w , equivalent to a Jordan block, has all the 
eigenvalues equal, it follows from the Perron- Frobenius theorem that the 
matrix must be reducible, more precisely applying recursively the Perron- 
Frobenius theorem to the square blocks of the equivalent triangular matrix, 
it is easy to show that the most general form for such a w is 



w 



P^UP = P^ 
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( X Ui2 U 
A ^23 

A 

V 



... \ 

A I 



P 



(24) 



For a Jordan matrix J, we also find that 



(fo-a) A^""" if a<h and b- a <n 
otherwise 



(25) 



and this implies that the susceptibility Xab can be written as (supposing for 
simplicity that the weight matrix can be put in a Jordan canonical form 
consisting of just one Jordan block) 



Xab = Rac E NiL)Jcd R'db ~ E RacR, 



_,T{-f + d-c) 



1 



db 



d>c 



{d-c)\ \<^' 



\\og{z. 



XK)) 



7+d— c 



(26) 



where we used the asymptotic form of N{L), the Stirling formula for the 

factorial, the formula (ESI) and Y,t x^L°'~^ ~ / . Now if we call A 

(losi) 

the biggest A = d — c for which I]c=i'~^ RacRc+A b fi^d the critical 

behaviour 



Xab 



T{i + A)E:z?-^RacR;l 



c+A,b 



1 



A! A'^ 



log{zXK)) 



7+A 



(27) 



From the formula (p^) and from w"' = P ^If^P, it follows that the Xaa have 
the usual critical exponent while the possible unusual critical exponents can 



and 



where both U and V are two upper triangulated invertible matrices with all the diagonal 
element equal 
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be find in Xah with a ^h. We can understand intuitively this phenomenon. 
Consider the oriented one-colour case (|l^, where = 0, then all the 
polymers beginning and ending with the same colour (±) are just made 
of monomers of the same colour as the ends. On the contrary a polymer 
beginning with + and ending by — is made gluing a + polymer and a — 
polymer, in such a way that the sum of the length of the two polymers is 
equal to the fixed total length, and this can be done in as many ways as it is 
the total length. 

In conclusion the critical behaviour of the non oriented coloured polymers 
is completely equivalent to the behaviour of the usual polymers, while the 
behaviour of the oriented (coloured) polymers can be different, even if there 
is just one set of critical exponents from which the others can be derived. 
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